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Abstract 

The NGT field equations with sources are expanded first about a fiat Minkowski 
background and then about a GR background to first-order in the antisymmetric part of 
the fundamental tensor, From the general, static spherically symmetric solution of the 
field equation in empty space, we establish that there are two conserved charges m and £^ 
corresponding to the two basic gauge invariances of NGT. There is no direct contribution 
to the flux of gravitational waves from the antisymmetric, gi/j,,^], sector in the linearized, 
lowest order of approximation, nor in the non-linear theory. It is demonstrated that the 
flux of gravitational waves is flnite in magnitude and positive deflnite for solutions of the 
fleld equations which satisfy the boundary condition of asymptotic flatness. 
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1. Introduction 

An analysis of the properties of gravitational radiation in the nonsymmetric grav- 
itational theory (NGT) is given, based on the field equations with sources. In previous 
work-*^, an analysis of the radiation problem in NGT, showed that the flux of gravitational 
radiation was the same as the leading order contribution in general relativity (OR). This 
result was shown to follow from the axisymmetric, time-dependent vacuum field equations 
in a first-order expansion about a curved GR background^. From an expansion of the 
exact, axisymmetric time-dependent field equations in inverse powers of r, the same result 
was obtained^. 

In the following, we shall analyze the problem using the NGT field equations with 
sources. In Sect. 2, we review the Lagrangian density and the field equations of NGT 
and in Sect. 3, we consider the results that follow from an expansion of the theory about 
Minkowski spacetime. Then in Sect. 4, we study the consequences that follow from the 
general, static spherically symmetric solution of the empty space field equations. We find 
that only two conserved charges occur in NGT, the mass m and £'^. These charges are 
conserved by virtue of the two basic gauge invariances of NGT, namely, diffeomorphism 
invariance and a invariance of the Lagrangian density. The other source tensor Tj^j^j 
contains only geometrical mass coupling. Indeed, when £^ = 0, NGT reduces to a purely 
geometrical theory of gravity involving only the gravitational charge m. It is found that in 
the flat space linearized theory, there is no direct contribution to gravitational radiation in 
the wave-zone. Only the GR quadrupole radiation manifests itself in this approximation. 

In Sect 5, it is shown that by using a spin-projection analysis, there are no propagating 
ghost modes. In Sect. 6, a calculation of the flux of gravitational radiation for plane waves 
is shown to be positive definite. 
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In Sect. 7, we expand the antisymmetric field equations about a classical curved GR 
background. We find that to first-order in the skew part of g^^, the fliix of energy is finite 
and positive for solutions of the field equations which satisfy the boundary condition of 
asymptotic flatness. As in the linear approximation, only the GR quadrupole radiation 
makes a contribution to the fiiix of gravitational waves. 

An important aspect of the results confirming the physical consistency of NGT, is 
that at no time is there any need for a gauge invariance in the skew g^i, sector of the 
theory. Indeed, such a general gauge invariance does not exist in the theory; the only 
invariances are general coordinate invariance and an Abelian gauge invariance, associated 
with the Lagrange multiplier field W^. This work demonstrates that NGT is a consistent 
theory of gravity, avoiding the criticisms of Damour, Deser and McCarthy^" 

2. NGT Field Equations with Sources 

The Lagrangian density with sources, in NGT, is given by^~^^: 

where g^''' — ^/—gg^'^ and Rij^y(W) is the NGT contracted curvature tensor: 

RAW) = - l{w^f3,u + wSf,,,) - w^^w^p + w^^^w-^, 

defined in terms of the unconstrained nonsymmetric connection: 

2 

where = W^^^^^ = 1{W^^ - W^^). This equation leads to: 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



The contravariant tensor g'^'^ is defined in terms of tlie equation: 

9'"'9a.=9''^9.a^Si^. (2.5) 
The NGT contracted curvature tensor can be written as 

R^4W) = R^4T) + ^W[^^,], (2.6) 

where i?^j^(r) is defined by 

The Lagrangian density for the matter sources is given by (G=c=l): 

Cm = -Stt^^^T^. + Y^^^S^. (2.8) 
Our field equations are given by 

G^AW) = SttT^,, (2.9) 

g[/^-] ^ = 47rS'^, (2.10) 



where 

a.... = /?..., - 
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G/^i, — R/j.^ — -z9nv'R- (2.11) 



The variation of the W connection gives 

2 Att 

+ g,f^w^, + g'^^w^;, - ^^^wp^ + -5:;g^^w^p + — (8-^5^ - s^5^) = o. (2.12) 

These equations can be written in the form: 

9^iu,a - 9pu^1j,a - 9t^pKu = 0, (2.13) 



where 

A^. = r^. + i^^.. (2.14) 

Here, D^^j^ depends only on and g^u and is defined by 

All 

9pvDP,a + ^MP^ai. = — i^SP{g^^gp^ - g^pgau + g^LuQiap])- (2.15) 

The variational principle yields for invariance under coordinate transformations the 
four Bianchi identities: 

[g'^^G',.(r) + g^'^G'.p(r)],, + g'-^pG.^iV) = 0. (2.16) 

The matter response equations are 

I (^apT'^" + 9paTn,a ' ^^a/3,pT"^ + ^^[p-^S^^ = 0. (2.17) 

From the invariance of JC under the U{1) gauge transformation: 

W'^ = W^ + X,^, (2.18) 

we obtain from Neother's theorem the identity: 

g[^'^],^,,^47rS^,^ = 0. (2.19) 

An important consideration, in determining the fundamental properties of NGT, is 
the possible invariances allowed in the theory. The rigorous NGT Lagrangian density 
contains the diffeomorphism invariance, associated with the invariance under general co- 
ordinate transformations, which leads to the existence of the four Bianchi identities (2.16). 
The only other invariance of the Lagrangian density is the Abelian invariance under the 
gauge transformation (2.18), from which follows the one identity (2.19). There is no rig- 
orous, general invariance associated with the g[^^], and there are no further identities 
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associated with such an invariance. However, as we shaU demonstrate, there is no need 
of any gauge invariance associated with g[^v']- This will prove to be true in the linear 
approximation to the theory, as well as in the higher-order nonlinear theory. 

3. Field Equations in the Linear Approximation 

Let us consider weak fields in NGT, whereby we expand ^^i^ about fiat Minkowski 
space^^: 

g,j,u ^ Vtxiy + h^i,, (3.1) 

where <^ 1 and 7]^,^ is the Minkowski metric: r/^i^ = (— 1,— 1,— 1,+1). We shall solve 
the field equations to lowest order in /i^j^. Raising and lowering is done using ry^j^, so that 
from (2.5), we obtain 

gt^i^ = r]'"' - hf"" + 0{h^), (3.2) 
h''^ = r^^'^^'^Kx- (3.3) 
Solving for A^^^ and D^^^ using (2.13) and (2.15), we get 

ri. = Iv^'^iK.,, + V,. - K,,a) - - (3.4) 

The full Lagrangian of the theory in the linear approximation is given by 

C = Cgr + Cs, (3.5) 

where 

JCgr = -Ih^^"^ □ - + i(/i('^«),« - lh'^)h^ + 87r/i(^.)T('^-), (3.6) 
and 



+ ^S^S^ - + 87r/i[^,]T['^-l. (3.7) 
Here, h = 'q^^h^^i, and □ = d^d^. The field equations to lowest order are: 

h^^pf =4.ttS^, (3.8) 

-^(n/i,^-/i(,,),^'"+/i,^,-^/i[^,],,''^ + ^/i[,,],^'") = ,]+87r(T^,-^r/^,T), (3.9) 

where T = v'^''T^„. 

By choosing the coordinate conditions: 

hia^y^^h^, (3.10) 

the field equations become (3.8) and 



nV.) = -167rT(^,), (3.11) 



= -^^If^M - Y'^t'^''^] + 167rT[^^], (3.12) 



where 



In the linear approximation, the skew Lagrangian Cs is not invariant under the gauge 
transformation: 

However, the presence of the Lagrange multiplier W^, in Cs, is of crucial importance, for 
it guarantees that unphysical ghost pole contributions will not occur in the theory. If the 
Lagrange multiplier were absent, then it would be vital to have invariance under the 
transformation (3.14) to avoid any unphysical propagating modes. 



It is important to distinguish between those degrees of freedom that propagate as 
physical dynamical modes, and those that do not couple dynamically to the matter sources. 
We shall therefore define longitudinal and transverse projection operators^: 



^ixu — ^ 1 Pjxv — Vnf (3.15) 



and make the general decomposition: 

h[ij,u] = tt/u,!/ - a^/./x + ^fj,uKxP^'^'^\ (3.16) 

and 

T[^.]-i^[MH+WA<^^"'^', (3.17) 

where and are a vector and a pseudovector, respectively. 
We choose the gauge condition: 

a^'^ = 0, (3.18) 

and find that 

hll] = P'^'^Pu^hap] = e^-A/^f'^'^l, (3.19) 

^S] = P^'^P.^^ap] = 0, (3.20) 

h^l] = P^rP^fho^P] = - (3-21) 
Using the gauge condition Wa'°' = 0, we have 

Moreover, we find that 

= e^,«A Jt'^'^^ (3.23) 



= 0, (3.24) 

It follows that 

The field equations with sources now become 



Ti'Z' =0. (3.26) 



/.g5j'^ = 47r5^, (3.27) 

□ /ij;^] =167re^,«A^t"'^^ (3.28) 

n^pM^] = - y%,.] + 167rX[^,.]. (3.29) 
By using (3.18), (3.21) and (3.27), we have 

na,^ = A7rS^. (3.30) 

Then, (3.21) and (3.30) give 

□ /ig:]=87r%,.], (3.31) 

and from (3.29), it follows that 

= -87r%,.] + 127rK[^,,]. (3.32) 

Thus, from (3.28) and (3.31), we see that the auxiliary (Lagrange multiplier) contribution 
^[/u,i/] can be eliminated from the field equations in the presence of sources. 

4. General Static Spherically Symmetric Solution 



In order to establish unequivocally the number of conserved "charges" in NGT, we 
shall study the consequences of the general, static spherically symmetric solution of the 
NGT field equations (2.9) and (2.10) in the absence of sources: S^^ = T^^, = 0. 
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In the case of a static spherically symmetric field^"^, Papapetrou has derived the 
canonical form of g^i,: 



/ -a{r) w{r)\ 

-Pir) f{r)smd 

-/(r)sin6» -/?(r)sin^^ 

\-w{r) 7(r) / 



(4.1) 



A general solution of the NGT field equations in vacuum was obtained by Vanstone^^, and 
is given by 



4(1 + 62)y 



^1 



(In y - a) 



Ay 



a 



1 



P + /3 



— y, w 



(4.2) 



where y is an arbitrary function of r, a is a complex constant, bi = 1 + is, and b, s, X and 
£^ are real constants. 

Consider now the asymptotically flat boundary conditions as r — > oo: 



CK — > 1, (3 r"^ , 7^1, w — > 0. 



(4.3) 



Vanstone proved that the only possibility to satisfy these conditions is for / to vanish 
identically. If we adopt the weaker condition that //r^ — > asr — > oo, and that the 
solution (4.2) must reduce to the Schwarzschild solution when gi^v'] = 0, then we get 



y = l 



2m 
r 



A = 4m^ 6 = 0, a = 0. 



(4.4) 



where the radial variable is only equivalent to the r variable in the GR Schwarzschild 
solution in the limit P ^ as r — > oo. 

We now put ^/h = /J^ + ii^, and obtain = ±[1/2 + l/2i/(l + s'^)]^/'^, v = ±s[2 + 
2y(l + s2)]-V2. Then we find that 

2m^[sinh^sinr/ — 6'(cosh^cos?7 — 1)] 



/ = 



(1 - — )(cosh^ - cos?7)2 



(4.5) 
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_ 2m [ssinh^sinr/ + cosh^cos?7 — 1] , , 

(1- 2^)(coshe-cosr/)2 ' ^ ■ 

«; = ±^, (4.9) 

where ^ = //Iny, 77 = zvlnj/. In the hmit that s = £^ = 0, we obtain the Schwarzschild 
solution: 

«=(i-v)"'' = ^ = i-v- 

Expanding (4.5)- (4.8) for small values of s, we get 

f=1^U-?^ + .X (4.11) 



3 V ^ 

. = .^(1-^....), (4..) 

The first important consequence of the general, static spherically symmetric solution 
is that there are only two conserved charges in NGT, the mass m and i'^. These are both 
conserved by virtue of the Bianchi identities (2.16) or (2.17), as in GR, and the NGT 
conservation law (2.19) with the identification: 



J d^xS^. (4.15) 



Thus, only two gauge invariances play a role in NGT, the diffeomorphism invariance and 
the invariance associated with the NGT Lagrangian density under the U{1) transformation 
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(2.18). There is no need for any extra gauge invariance associated with the g[ij,u] sector, 
for there does not exist a conservation law of a new charge in the theory besides m and 

The second significant consequence of the static solution is that in the linear approx- 
imation, neglecting contributions of order h'^^j^^y we see from (4.11) that h[23] is zero to 
lowest order. From (3.28), it follows that to lowest order of approximation (neglecting 
terms quadratic and higher order in m) in the static limit: ~ ^ therefore in 

general: /i^^j,] = 0. Then, to lowest order we have 

J'' = 0, %.]=i^[^,.]. (4.16) 
It now follows from (3.32) that to lowest order of approximation: 

in the wave-zone. Because T^^^^ is a compact, localized source tensor, then from (3.17) and 
(4.16), it follows that has compact support and vanishes in the wave-zone together 

with 5''^ and (4.16) follows immediately. 

In the wave-zone at infinity, the skew vacuum field equations are given by: 

ht'^/ = Q^ (4.18) 

n^fH=0- (4-19) 

The solution for /t(^i/) coincides to lowest order with the GR solution^^: 

/.(,.)(x,t) = -4 / dV ^(^-^^"^-'_''"~"^'^ (4.20) 
J 1-^ I 

To this order /i^^j = and the solution of equation (3.31) is given by 

fta(x,t) = 2 / A' ^""-''t-"'r"'" - (^-21) 
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5. Spin Projection Analysis of the Linearized Theory 

Let us write the skew part of the Lagrangian Cs in the form^^: 

>Cs = ^ ^ (f)AOAB(f>B, (5.1) 
A,B 

where (pA = (^[/ii/], W/j) and Oab is the wave operator. Using the methods of ref (16), we 
can decompose the fields into subspaces with spin-parity and invert the operator Oab 
to obtain the saturated propagator: 

n = - E SaO^'bSb, (5.2) 

^A,(I>B 

where Sa = {L[i^„],S,^) and 

L[^l,y] = T[^u] + %,!/]• (5-3) 
Expanding the operator Oab yields 

^s= E (5.4) 

ll^A,(l>B,i,j,J^ 

and 

n = - E a-/''^{JnSAP;!'/{JnABSB. (5.5) 
where af^^{J^) are matrix coefficients. A calculation of the spin-projection operators yields 

^^■(1+) = li0^.a0.p - e^J^p), (5.6) 

- \ _ ( \ beta — ^ixP^va — ^va'^ij.p + ^vP^ixa) -^{ki^^fia ~ k^dj^a) \ 

Pij{Q-^)^uj^p, (5.8) 

where 

— Vi^a ~ k^kak , w^q, = k^kak , k^ — k^ik ) ^ . (5-9) 
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These operators are ortho normal and complete within W^). The Lagrangian 

Cs now yields the matrix coefficients: 

a(l+) = Ik', (5.10) 



= ( ^(;,2)i/2 ' ) ' (5-11) 



a(0+) = 0. (5.12) 

The result (4.12) follows from the gauge invariance under the transformation (2.18) and 
the conservation equation (2.19). Inverting (4.10) and (4.11) yields 

a-\l-)=( ' -f.i^')'^']!. (5 14) 
A calculation determines the saturated propagator to be of the form: 

n = -pT['^^]p(l+)^,«^T["^] + —S"{k,9^a - k^9.a)L^'"'^ 

q 

;L^^''\k^d,a - k.e^o)S^ + -(5«^«/3^^). (5.15) 



Here, we have used the condition k°'9ai3 = to remove terms of the form P{l~^)ij,v>ai3 

We see that the massless spin 1~ ghost particles do not propagate, and only contact 
terms occur in this sector, whereas the spin 1+ sector is neither a ghost particle nor a 
tachyon. For the rea/ version of NGT, we see that the theory is completely ghost free, even 
with a nonvanishing T^^u] ■ 

The nonpropagating spin 1~ sector corresponds to the /i^^j components, while the 
propagating spin 1+ sector corresponds to the transverse hj'^-^ components of In 
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Sect. 6, we will prove that the ^^"^j do not propagate in the wave-zone. Although the 
skew field equations are not invariant under the transformation (3.14), there does exist a 
restricted gauge invariance under the transformations: 

ShfiL] = ^M,- - ^hj;^.] = WaC''^'^^ (5.16) 

provided that and satisfy: 

□e^ = 0, e^'^ = 0, DCm^O, Cm''' = 0- (5.17) 

This restricted gauge invariance reduces the three degrees of freedom of the propagating 
spin 1+ sector to one spin 0+ degree of freedom. It should be stressed at this point that 
we do not need to use this restricted gauge invariance to remove potential ghost modes in 
the linear approximation, for such ghost modes simply do not exist. 

A similar result occurs for a spin 1~ field with a scalar Lagrange multiplier 
(/), if we replace with and with ^ . The spin 0"*" ghost particle becomes 

nonpropagating because of the Lagrange multiplier cf). 

If we had included a mass term in £5 of the form cigf'^^^5r[^i/], then we would be forced 
to have an unphysical particle spectrum, which contains either ghosts {d < 0) or tachyons 
(d> 0). Only for the massless theory are we guaranteed the absence of ghosts, due to the 
presence of the Lagrange multiplier W^. The suggestion by Damour, Deser and McCarthy 
that a mass term be introduced into NGT to avoid bad asymptotic behavior would 
necessarily lead to an unphysical theory. 
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6. Plane Wave Solutions in the Linear Approximation 

The plane wave solutions for are identical to those in GR. The solution of Eqs. 
(4.19) is 

hfiZ] = e£X]exp(i/cAa;^) + el^Jjexp{-ikxx^), (6.1) 

where /c^/c^ = 0, 

k^e[^.] = 0, (6.2) 

and to this order hj^-^ = 0. 

Let us now calculate the energy-momentuum tensor of the plane waves. This is done 
using the energy- momentum pseudotensor t^i, of NGT. We have to lowest order^^: 

V = ^«nr)-^r7^.i?^')(r)), (6.3) 

where i?p^J(r) is i?^iy(r) to 0{h?). We only need the average expression, < >, since 
this is what is measured. Averaging over space and time in a region much larger than 
|k| ^, we have 

This expression is positive definite for the real (or hyperbolic complex^^) version of NGT. 
Here, we have taken into account that is zero to lowest order of approximation in 

the wave-zone. 

7. Generation of Gravitational Waves in the Linear Approximation 

The analysis of the gravitational waves generated by the the GR source, in 
the linear approximation is well understood^^. Let us consider the possible existence of 
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radiation generated by the NGT sources, and T^^u] , in terms of their Fourier components 
(x, t) and Tj^j^] (x, t) in the hnear approximation: 

/•OO 

5^(x,t)= / dwS^'{x,oj)exp{iojt) + CC, (7.1) 
io 

/•OO 

%^](x,t)= / da;K[^,^](x,a;)exp(za;t) + CC, (7.2) 

where CC means complex conjugate. 
In the wave-zone, we get 

|x-x'| =r-x-x' + ... (7.3) 
We have /i^^j = and from (3.31), we obtain 







|x-x'|)] 




X — x' 





-a.^ydv ------ J. (7.4) 

By using the expansion (7.3), we get 

h^l] = ef^l] (x, u;)exp{ikxx^) + CC, (7.5) 

where 

efi = l[KS^O^,u;)-k^S.ik,oj)]+Oil/r^). (7.6) 
When we substitute (7.6) into (6.4), we find that 

< V) >= ^K^'^^'^w - ^l''!^)]- (7.7) 

Thus, the e^^j do not contribute to the gravitational wave flux through the source -Sj^^j^], 
and this demonstrates that /j.^^^ does not propagate. This agrees with the conclusion, 
reached in Sect. 4, that the spin 1~ sector does not propagate. Since the NGT pseudovector 
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source J'^ = to linear order, the gravitational wave flux is the same to leading order as 
in GR. This coincides with the results obtained previously-*^"^. We observe that < > 
is positive definite, and the linear approximation is physically consistent. 

8. NGT on a Curved GR Background 

We shall expand the NGT field equations in powers of taking the symmetric 

part g{fj,iy) to be an exact GR background, g'^^. To lowest order in 

The field equations with sources become in lowest order: 

4 1 Stt 

D'^Flxva + "^RoR^J^vhlaP] = -^^^[m,^'] + 167rT[^j,], (8.3) 

= 47r5^, (8.4) 
where D"' is the GR background covariant derivative and 

F^up = h[^y]^p + h[^p]^^^ + (8.5) 

We observe that these equations are not invariant under the gauge transformation (3.14), 
as is expected, for this invariance does not occur in the theory. 

Because of (3.17) and (8.5), we can rewrite Eq. (8.3) in the form: 

D'^Dah^^,] + ARGR''^'^,h[af3] = W[^,u] + IGire^.a^J^'^ (8.6) 

where 

Wl^.M = - V[M,.] - ^S^f^M + 167rK[^,.]. (8.7) 
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By taking the curl of (8.6), we obtain the field equation: 

D'^DotF^j^yp + -D{p(5[^z^]} = 0, (8.8) 
where {...} denotes the curl. Moreover, we have 

= 4i^Gi^"/>[a/3] - 167re^,«;3 J«'^. (8.9) 
On the other hand, we could take the divergence of Eq. (8.6) and get: 

D'^D^W,, = -D^{8RGR%^h[ap]), (8.10) 
where we have used the gauge condition: 

D^Wf, = 0. (8.11) 

Eq. (8.10) represents four integrability conditions that must be satisfied, once the 
solutions for and /ij^j,] are obtained. The boundary conditions are imposed on 
and h[^jy-^ by demanding that g^^ r]^^, and 0. 

It has been proved by explicitly solving the system of six equations: 

D''h[^,] = 47tS^, (8.12) 

D°'DaFfj^j,p = (8.13) 

that consistent solutions can be obtained for the six which satisfy the physical bound- 
ary condition of asymptotic fiatness at infinity^"'^. It is important to observe that any h^^u] 
solutions of (8.12) and (8.13) must satisfy the primary equations (8.6). The same is true 
for any solution of (8.10)^-2°. 
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The skew contribution to the stress-energy pseudotensor, t^^^'^\ can be expressed in 
the form^: 

The GR Riemann tensor must have dimensions of [length]"^ and contains at most two 
time derivatives. The dimensionless quantity, M, (where M is the mass function) appears, 
together with the quantity M, and therefore the Riemann tensor behaves in the wave-zone 
as ~ 1/r. The rate of energy loss for an isolated body, in NGT, is given by 

where the integration is over a sphere of radius R in the wave-zone at infinity and Ui is 
an outward pointing unit vector. At infinity, the first term in (8.14) is positive, while 
the second term does not contribute, since in the flat spacetime limit is zero to 

lowest order, as was shown in Sect. 4. The third term goes to zero at least as ~ and 
does not contribute to the energy-momentum fliix, while the fourth term has the form: 
dadpN^f^°'^^'^f^\ making no contribution. Therefore, the total energy-momentum flux in the 
wave zone is finite and positive definite. 

In the case of an axisymmetric, time-dependent solution of the NGT field equations 
to first-order in on a curved GR background^, and for the generalized Bondi, van 
der Burg, Metzner and Sachs (BBMS) solution^'^-*^, the calculations show that only the 
quadrupole GR radiation manifests itself in the energy-momentum fiiix, for solutions that 
obey the asymptotic flatness boundary condition at future and past null infinity. This 
agrees with the result obtained in the linear approximation in fiat spacetime. However, as 
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in the case of GR, it is not possible in the generahzed BBMS solution to interpret directly 
the nature of the near-zone NGT sources. The properties of these sources can only be 
inferred from multipole moment expansions and their effects on wave solutions at future 
and past null infinity. It should be stressed that the NGT source in the axisymmetric, 
reflexion symmetric solution represents a generic source in NGT, for it possesses all the 
GR and NGT moments of a general gravitational wave source. 

9. Conclusions 

The proof of the physical consistency of NGT does not depend on the existence of 
any form of gauge invariance in the g[p,,y] sector, for there exist only two conserved charges 
consistent with the two gauge invariances, namely, diffeomorphism invariance and the 
U{1) gauge invariance that guarantees the conservation of the charge For an isolated 
body, the physical boundary condition demands that in orthonormal, Cartesian coordinates 
g^j^ — > ?7pi^ and g^^^-^ — > as r — > oo, which in turn requires that for the static spherically 
symmetric NGT solution in vacuum, ^[23] = to the lowest order of approximation. In 
the higher non-linear orders of approximation there is no direct contribution to the flux 
of gravitational waves from the skew sector. It then follows that NGT gives physically 
sensible answers for the radiation of gravitational waves. 

Axisymmetric time-dependent solutions have been obtained, which show that the 
expected consistency of the field equations does hold. A solution to the axisymmetric case 
was derived to first-order in in an expansion about a curved GR background^. The 

cross-coupling term between the background GR Riemann curvature tensor and /ij^j^j does 
not cause any inconsistency, provided the field equations are correctly solved using the 
boundary conditions imposed on g^^,, namely, that the solutions are asymptotically flat 
at past and future null infinity^. The exact axisymmetric NGT vacuum field equations 
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were expanded in inverse powers of r, and it was found that the gravitational wave Aux 
of energy was positive definite^. This result extended the proof of the positivity of the 
gravitational wave flux in GR, obtained by Bondi et al. and Sachs^-*^'^^. 

Because there are only two conserved charges, m and in NGT, it can be proved 
that there cannot exist any direct contribution from Qi^i,] , which has the form of a retarded 
or advanced wave with the asymptotic behavior g^^^^-^ pa 1/r. This follows from the fact 
that has the dimensions of a [length]^, whereas m (and the electric charge e in electro- 
magnetism) has the dimensions of a length and can generate a 1/r retarded wave solution 
and quadrupole radiation^"^. 

When the NGT charge £^ = 0, the static spherically symmetric solution only contains 
the mass m and the theory becomes purely geometrical. Test particles will follow geodesic 
equations determined by the Christoffel symbols of NGT^^. 

From Eq. (8.2), we see that to first-order in h^^j^j, a positive energy theorem follows 
trivially using the methods of ref. 25. However, further work must be carried out to derive 
a general, rigorous positive energy theorem in NGT. 
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